Monochromatic path crossing exponents and graph connectivity in 2D percolation 
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We consider the fractal dimensions dk of the fc-connected part of percolation clusters in two 
dimensions, generalizing the cluster {k = 1) and backbone (k = 2) dimensions. The codimensions 
ifc = 2 — dfc describe the asymptotic decay of the probabilities P{r, R) ~ (r/RY*' that an annulus of 
radii r and i? ^ 1 is traversed by k disjoint paths, all living on the percolation clusters. Using 
a transfer matrix approach, we obtain numerical results for Xk, k < 6. They are well fitted by the 
Ansatz Xk = ^1^^ + + {1 - k)C , with C = 0.0181 ± 0.0006. 
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Percolation is a classical model of statistical mechanics 
and plays an important role in the study of disor- 
dered systems It is also one of the simplest models 
displaying a critical point. In two dimensions, exact val- 
ues for a variety of critical exponents have been found 
over the last two decades, and quite recently many of 
them have been confirmed by rigorous probabilistic ar- 
guments Q. Most of these exponents can be defined 
through the fractal dimensions of suitably defined sets at 
the percolation threshold. 

In the present Letter we shall be concerned with an 
infinite family of critical exponents whose exact values 
remain unknown to this date. These exponents charac- 
terize the connectivity structure of the percolating clus- 
ter (s) at criticality. 

Following Tutte we define a graph to be k- 
connected (for fc > 1) if no separation into disconnected 
subgraphs is possible by eliminating at most fc — 1 vertices 
along with their ingoing edges. It is easy to see that we 
may equivalently require any two vertices in the graph to 
be connected by (at least) fc disjoint paths. 1-connected 
graphs are simply the percolation clusters, and to inquire 
into the connectivity structure of a given cluster, we may 
decompose it into its largest 2-connected components ||] 
(better known as 2-blocks, or "blobs" , in the percolation 
literature), and so on. Tutte has shown that the decom- 
position of a 2-connected graph into its largest 3-blocks is 
unique g. 3-blocks are relevant for applying Kirchhoff's 
laws to resistor networks, and are useful for analyzing the 
performance of certain algorithms . 

To better study the transport properties of percolation 
clusters, we henceforth consider critical percolation in a 
large square of linear size L, and we specialize to clus- 
ters that connect to the boundary of the system. In the 
limit L — > oo, the boundary becomes the "point at infin- 
ity" , and the above definition states that a given vertex 
is fc-connected if it is connected to infinity by (at least) 
fc disjoint paths within a percolating cluster. Following 
Ref. , we shall call the set of fc-connected vertices the fc- 



bone. The 2-bone is of course nothing but the (geometri- 
cal) backbone, i.e. the part of the percolation cluster that 
sustains a non-zero current, when a voltage difference is 
applied between its two terminal points. (As usual we 
disregard rare Wheatstone's bridge-like arrangements.) 

We shall here be interested in the fractal dimension 
of the fc-bone, assuming its mass to change with system 
size like L'^'' . The cluster dimension rfi = f| has been 
known exactly for a long time [0-^. The backbone di- 
mension has recently been related to the solution of a 
partial differential equation which however appears 
to be intractable, even numerically. Still, numerical es- 



timates are available from Monte Carlo |10| and transfer 
matrix methods |l|]: da = 1.6431 ± 0.0006. After the 
completion of this work, a first estimate for appeared: 
da = 1.2±0.1 1^. Actually, this result was obtained 
from block-decomposition of clusters and backbones, but 
for reasons of universality we expect it to apply to the 
3-bone as well. Also, the above definitions are stated for 
site percolation, but the exponents should be the same 
for bond percolation, with the clusters being separated 
by cutting edges rather than vertices. 

A useful alternative formulation of the fc-bone problem 
is obtained by passing to an annular geometry, limited by 
two concentric circles of radii r <C 1 and i? 3> 1 , by means 
of a conformal mapping. (This is permissible since per- 
colation has been proved to be conformally invariant [Q.) 
Interpreting the inner circle as the point which is a po- 
tential element of the fc-bone, and the outer as the point 
at infinity, we see that a given percolating configuration 
in the annulus contributes to the fc-bone if and only if the 
two circles are connected by fc disjoint paths on the per- 
colating cluster(s); see Fig. |l|. The fractal dimension dk 
of the fc-bone is linked to the scaling of the path-crossing 
probability Pfc(r, R) ~ (r/i?)^*" through the codimension 

= 2 - dk iTai. 
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FIG. 1. Annular geometry endowed with critical percola- 
tion (here shown in the continuum limit). The existence of 
two disjoint traversing paths on the clusters implies that this 
configuration contributes to the 2-bone. 

A more general class of path-crossing exponents can be 
defined from traversing configurations where some of the 
paths live on the percolating clusters (black paths), and 
the rest on the dual clusters (vi^hite paths). Interestingly, 
the corresponding critical exponents Xk only depends on 
the number of paths, fc > 2, and not on their colors, 
provided that both are represented jlj]. Their values 

x,^^{k'~l) (1) 

are known rigorously and differ from those of the 
monochomatic exponents Xk defined above. 

Some information on the Xk is provided by the inequal- 
ities 

Xk < Xk+i- (2) 

This inequality is valid since a configuration contribut- 
ing to Xk+i can be taken to have k black paths and one 
white path. Clearly, it then also contributes to Xk- Fur- 
thermore, as fc — !■ oo, the effect of a single extra path 
should be small and we expect it to be of the same order 
if one changes the color of one of the k existing paths; 
thus, Xk+i - Xk ~ Xk+i - Xk = 0{k). 

In view of Eq. (Q) we find the asymptotic result Xk — 
j^k^ + 0{k) as fc — > oo. In analogy with Eq. (^ it 
thus seems natural to conjecture that the spectrum Xk 
is quadratic in fc. From ii = ^ we then get 

ik = ^k^ + ^k+(l-k)C, (3) 
with C — 0.0181 ± 0.0006 from the numerical result on 

To check the conjectured form of the spectrum, Eq. (Q), 
we now turn to our numerical results. But first we must 
briefly describe our transfer matrix algorithm; it is a nat- 
ural generalization of the one used in |ll|. 

First we consider the annulus of Fig. Ill as a cylinder 
with circular space and radial time. We have tried several 
choices of lattices to discretize it. For practical applica- 
tions, it turns out to be best (see 0) to use a square 
lattice with a "light-cone" orientation, that is such that 
the periodic direction forms a 45 degrees' angle with the 



two axes of the lattice. We then define the discrete time 
slices such that they intersect the lattice at vertices only: 
we call L the number of such vertices (in units of the 
lattice spacing the period is then Ly/2). 

Next we define the basis of states on which our transfer 
matrix acts. A basis state is a collection of path configu- 
rations; a path configuration is the data of the positions 
of our fc paths at a given time, with possible additional 
"arches" to allow backtracking of paths, see Fig.|. Note 
that the encoding of states can be easily implemented 
as follows: basis states are encoded as sorted lists of 
path configuations, and path configurations are repre- 
sented by words of length L made out of the four letters 
{opening^ closing, path^ empty} and which contain fc let- 
ters path. The letters opening and closing define the 
backtracking arches. 

Ill" irr ir*i nr rn 
r*n 'iir 'in tii "iii 
"in i^n iri iir-^iir- 

FIG. 2. Path configurations for fc = 3, L = 5. 

The transfer matrix itself acts on a basis states by 
"evolving" all the path configurations it contains with 
a single configuration of bonds (i.e. percolating/non- 
percolating state of each bond) and recombining the re- 
sult into a single state, then summing over all configu- 
rations of bonds. Evolving a path configuration means 
considering all possible continuations of the paths from 
one time slice to the next, including possible appearances 
of new arches, or existing arches connecting to paths or 
other arches. We use sparse matrix factorization tech- 
niques to build up the entire transfer matrix; dihedral 
symmetries are quotiented over once a time slice has been 
completed. 

As in [pl| , it is convenient to allow "perpendicular tan- 
gencies" for paths even though they should be in princi- 
ple excluded; that is to allow two paths to touch at one 
vertex in the configuration whenever the tangent of ei- 
ther path is perpendicular to the transfer direction (we 
still exclude "parallel tangencies"). It is expected, and 
we have numerically verified, that inclusion of either or 
both types of tangencies and/or changing the lattice ori- 
entation does not alter the first finite-size correction of 
the eigenvalue of the transfer matrix, which yields Xk- 
However, the particular choice of allowing perpendicu- 
lar tangencies on the light-cone oriented lattice has the 
advantage of greatly decreasing the number of states gen- 
erated (see below) and improving the convergence prop- 
erties of the finite-size data. 

Finally, the following procedure is used to compute the 
matrix elements of the transfer matrix: Starting with an 
arbitrary basis state (e.g. the one consisting of the single 
path configuration path^ empty^~^) , one acts on it with 
the transfer matrix, stores the corresponding matrix ele- 
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merits, then considers all new basis states generated and 
iterates the procedure until no further new states are gen- 
erated . What we build this way is a submatrix of the 
transfer matrix corresponding to a stable subspace; this 
submatrix is in fact much smaller than the full transfer 
matrix, which is essential for practical applications. One 
then extracts its largest eigenvalue Afc(L), which yields 
the codimension Xk via the formula 



2^A,(i) = l-^+o(L-i) 



(4) 



We show in Tab. | the results for 2 < A: < 6. Since the 
memory and time requirements presumably grow facto- 
rially with L, we cannot push the calculation very far in 
L. It is however sufficient to estimate Xk, which is also 
given in Tab. | together with approximate error bars. 
The data for fc = 2 are in fact taken from jl^ , where the 
state space was reduced by exploiting that the case of 
two paths can be treated as an extra backtracking arch. 

For A: > 2, all the exponents Xk are consistent with, 
but less precise than, the conjectured spectrum ^ with 
C given by X2 ■ We also note that the inequality Xk < Xk 
seems to be satisfied, and if would be interesting if one 
could prove this. Another open question is the possi- 
ble rationality of the Xk] in this respect, the failure of 
computing these exponents by conformal field theory is 
particularly intriguing. 

For fc > 5, the dimensions dk — 2 — Xk are negative. 
Physically this means that fc-bones with fc > 5 become 
increasingly rare as the system size increases. Of course, 
the lattice model does not support a fc-block when fc ex- 
ceeds the coordination number of the lattice. However, 
the exponents Xk characterize the continuum limit, and 
are thus believed to be independent of the microscopic 
details; an appropriate lattice definition of the fc-bone 
for high fc is obtained by demanding that fc independent 
paths connect any small neighborhood to the point at 
infinity. 

Finally, we remark that the fc-bone problem extends 
to the Kasteleyn-Fortuin representation [Q of the q- 
statc Potts model (bond percolation being the limit 



g — > 1). Our numerical algorithm can straightforwardly 
be adapted to this case as well. 
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TABLE L Eigenvalues \k{L) of the transfer matrix and estimate of Xk for 2 < fe < 6. 



_k 2 3 4 5 6 

T 

4 0.718747415570 0.413598206498 0.121093750000 

5 0.775012703547 0.526618869796 0.257122218539 0.061523437500 

6 0.812529692986 0.603476157424 0.362299981029 0.153371684616 0.031005859375 

7 0.839330907375 0.658986646726 0.443031423565 0.237989873966 0.088905009155 

8 0.859432882632 0.700919030179 0.506272495802 0.310651059489 0.150977532764 

9 0.875067710677 0.556925756584 0.372225212770 0.210050339522 

10 0.263993624780 

ife 0.3569 ± 0.0006 0.77 ±0.02 1.33 ±0.03 2.1 ±0.2 3.0 ± 0.3 
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